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Abstract 

Low energy effective actions arising from string theory typically contain many scalar 
fields, some with a very complicated potential and others with no potential at all. The 
evolution of these scalars is of great interest. Their late time values have a direct impact 
on low energy observables, while their early universe dynamics can potentially source 
inflation or adversely affect big bang nucleosynthesis. Recently, classical and quantum 
methods for fixing the values of these scalars have been introduced. The purpose of 
this work is to explore moduli dynamics in light of these stabilization mechanisms. In 
particular, we explore a truncated low energy effective action that models the neigh- 
borhood of special points (or more generally loci) in moduli space, such as conifold 
points, where extra massless degrees of freedom arise. We find that the dynamics has 
a surprisingly rich structure — including the appearance of chaos — and we find a 
viable mechanism for trapping some of the moduli. 
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1 Introduction 



One important obstacle to the construction of realistic models of string phenomenology is 
the generic existence of moduli in the low energy theory [H El El HI El E]. The compact 
manifolds satisfying the string equations of motion generally come in continuous families 
whose parameters (controlling deformations of the metric and p-form fields) become scalar 
fields with exactly flat potentials in the 4 dimensional effective theory. If on-shell branes 
are present, there are also moduli that parametrize their relative positions and orientations. 
Additionally, irrespective of the compact manifold, low energy string models always contain 
the massless dilaton. 

These fields can cause a number of tensions with observation. For example, their presence 
in the early universe can modify the abundances of hydrogen and helium, a well confirmed 
prediction of big bang nucleosynthesis [7]. If there are exactly massless fields at late times, 
thermal and quantum fluctuations can cause time variation of standard model parameters 
and 5th force type violations of the equivalence principle [H [H [10]. On the other hand, 
fields with relatively fiat potentials can be useful in building models of slow-roll inflation. 
For these reasons it is important to explore the variety of moduli space dynamics as fully as 
possible. 

Proposals to mitigate the phenomenological issues mentioned above typically involve ei- 
ther the introduction of a potential for the moduli fields [121 [131 [IH [ISl [HI [IZ] , or suppression 
of their couplings to matter [TSl [T^l [20] • An example of the former strategy that has gener- 
ated much interest recently is to imagine the universe settling into a vacuum with nontrivial 
flux threading cycles in the compact extra dimensions [HI [12] . This induces a potential on 
the moduli space [121 [13] with a great number of apparently consistent vacuum solutions 
[211 [22] : perhaps an infinite number [23], or perhaps merely 10~^°° |24j^ . In light of statisti- 
cal analyses of this 'landscape' that find large numbers (but small compared with the total 
number of vacua) with small positive cosmological constant, some have reconsidered the 
anthropic framework to explain details of our particular universe fIT\ [2B] that have resisted 
more traditional approaches. 

Moduli space dynamics enter the discussion in two qualitatively different ways. Some 
approaches are aimed at avoiding anthropic reasoning through dynamical selection, while 

^It is worth keeping in mind that the landscape is not an estabhshed feature of string theory. A large 
number of apparently consistent vacuum solutions have been constructed, but their consistency has not been 
demonstrated beyond doubt, and there remain worries about whether they should be regarded as separate 
vacua of a single theory. ^2hi 
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others seek to bolster the anthropic approach by providing an underlying mechanism to 
populate a vast arena of possible universes. An example of the former strategy is the attempts 
to find a wavefunction on the moduli space [271 [2H] , whose square might be interpreted as 
a dynamical weight on the landscape. An example of the latter is the application in [21] 
of results on semiclassical tunneling [29] to the case of transitions between flux vacua with 
different cosmological constants. 

It was suggested in [301 E] that instead the problem lies with expanding about a single 
minimum of the flux potential and that the actual ground state is a relaxed superposition 
of all of the degenerate connected minima. 

The approach we take here is to study moduli dynamics directly in the low energy effective 
fleld theory, focusing not on generic loci in moduli space, but on neighborhoods of points 
where extra massless degrees of freedom appear. Section [2] motivates our choice of action by 
considering as an example the effective fleld theory resulting from a string compactiflcation on 
a Calabi-Yau 3 fold near a conifold point. Sections [3] and S] study two proposed mechanisms 
for stabilizing the moduli near these so-called extra species points, or ESPs. The idea of 
classical trapping has been considered in a number of contexts; early work in this direction 
includes [32] while some more recent examples are gases of wrapped strings and branes 
(see [53] and references therein), black hole attr actors |i34j, D-brane systems [35j, M-theory 
matrix models [36J, flop transitions [37] and conifold transitions [5Sl EZl [M] studied in a 
cosmological background. In the case of quantum trapping, previous works include [HI HJl 
H3l HU Hi]. In [m US] trapping is used to study trapped inflation and in the latter trapped 
quintessence. Here, in the conifold context, we go further than previous studies by taking 
account of both classical and quantum considerations, and of prime importance, establishing 
that the moduli systems we consider exhibit chaotic dynamics. Chaotic moduli evolution is 
something that had previously been suggested for moduli dynamics in flat spacetime; we use 
the Poincare 'surface of section' technique to put this idea on a flrm footing. Surprisingly, it 
turns out that an understanding of the chaos that appears in the Minkowski case is crucial 
to interpreting the dynamics when Hubble friction from the gravitational background is 
included. Speciflcally, we flnd the role of Hubble friction in classical moduli trapping to 
be less effective than previously believed, and that quantum particle production traps the 
moduli (or at least stabilizes their expectation values) far more efficiently. 

To facilitate a direct comparison with refs [Ml EH], we work with a 5d anisotropic back- 
ground. Our qualitative conclusions remain the same in the 4d isotropic case. 
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2 String Theory Motivation 



The low energy effective action we use in the remainder of this paper is relevant to many 
moduli systems. In this section we outline one context in which it arises — a Calabi-Yau 
compactification of M-theory|^ Readers primarily interested in the dynamics of the low 
energy theory, can skip to the result, which is equation (fTOl) . 

Our starting point is the low energy description of M-theory, which is the effective field 
theory of lid supergravity with action 



S 
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where F^ is the four- form flux F4 = dA^, and = SttGu = SttI^, with Gu the eleven- 
dimensional Newton's constant and Ip the Planck length. We want to consider compactifi- 
cation to 5D on a Calabi-Yau (CY) 3-fold /C with Ricci-fiat metric gab- Following [38j, we 
will take the lid metric to be of the form. 



ds' = -de + aitfdx' + h(tfdy' + gab{z)dz''dz^ 



(2) 



where the spacetime coordinates are = {t, Ob J Ob ^ Ob ^, y) and the CY coordinates are denoted 
by z"". 

Performing a Kaluza-Klein reduction to 5 dimensions and keeping only the massless 
modes results in an = 2 abelian gauge theorjjf containing; (see e.g. |47j ) 

1 gravity multiplet, 

(ea^,V'M/,^M)' (^ = 1,2) (3) 

hi^i — 1 vector multiplets. 



(2 = l,...,(/l(l,l)-l)) 



(4) 



and h2.i + 1 neutral hypermultiplets. 



{C\A\) 



(Z = l,...,2(/i(2,l) + l)) 



(5) 



^The full derivation contains many technical details of little relevance here, but can be found in [38l[37l[39] . 

•^Here we are using the language of 4D supersymmetry. In 5 dimensions, iV = 2 is the minimal nontrivial 
supersymmetry, and so is sometimes referred to as A'^ = 1. The unambiguous statement is that there are 8 
supercharges. 
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where hij = dim Hl^-' {JC) are the Hodge numbers of JC. The scalar fields in these multiplets 
parametrize the Kahler and complex structure moduli spaces of /C controlling its size and 
shape. Since the gauge theory is abelian and none of the hypermultiplets are charged, the 
potential for the scalars is exactly fiat, which is just to say that they are moduli. 

The above is not the full story, however, as the metric on the hypermultiplet moduli 
space is geodesically incomplete]^ i.e. a classical scalar field can reach the boundary in finite 
proper time. Originally this was seen as a proble nfl because boundary points correspond 
to singular configurations of JC that arise when cycles on /C shrink to zero volume. This 
is potentially troubling because the geometrical singularity gives rise to various divergent 
coefficients in the effective action. Type II string theory and M theory have a remarkable 
way of resolving the issue [H] [19] . Briefly, the singularity in the action can be interpreted as 
arising from incorrectly integrating out brane wrapping modes that are usually massive (since 
their mass is proportional to the volume of the cycle they wrap), but become massless at 
these special points in moduli space where cycles collapse to zero volume. Calculations show 
that by including these new massless modes in the effective action, the physical singularity 
is resolved. It is then natural to seek the explicit form for the low energy effective action 
near such a geometrical singularity [38j. The crucial point for our purposes is that the brane 
states are charged under some of the vector multiplets, so the modiflcation required to the 
4D theory is to introduce some extra charged hypennu.t.p.etsB It then follows from the 
extended supersymmetry that several of the scalars acquire potentials. The reason is that 
the scalars in the charged hypermultiplet must couple to the scalars in the vector multiplet 
corresponding to the charge they possess. 

Unfortunately, little is known about the fully quantum corrected hypermultiplet metric^ 
so it is not yet possible to write explicitly the 5D action corresponding to a compactiflcation 
of type II string theory on a given Calabi-Yau 3-fold JC. The best we can do at present is to 
note that N = 2 supersymmetry implies that the hypermultiplet metric has holonomy group 
precisely Spiriu) ■ Sp{l). Spaces admitting such metrics are termed quaternion-KdJiler, 
a somewhat confusing name since as manifolds they are quaternionic, but not in general 
Kahler, or even complex. Locally symmetric spaces (i.e. Lie group quotients G/H) admitting 

^The same is true of the vector multiplet moduli space, but we won't discuss that case here. 
^See for example the introduction of ^35]. 

^Of course the dimension of the parameter space of the scalars will increase, so the hypermultiplet 'moduli 
space' metric will change as well. This change is necessarily nontrivial, since simply taking a product of two 
quaternion-Kahlcr metrics does not result in another quaternion-Kahler metric. 

^An indication of the current status of quantum corrections to hypermultiplet metrics can be found in 

m 



5 



quaternion-Kahler metrics were studied by Wolf [5T] , who classified them into a few (infinite) 
families. For these examples one can find analytic expressions for the metric, but it is 
unknown whether any of them correspond to the moduli space that would arise from a 
genuine Calabi-Yau compactification. 

We will now outline the derivation of the full 5d effective theory including the additional 
massless degrees of freedom. A more detailed discussion of the construction can be found in 



[38] and [37]. There it was shown that one can pick the Wolf space u(^n\j)xu{2) ' truncate 
the action so that only two scalar fields are nonzero; one from a vector multiplet, and one 
from a charged hypermultiplet. The result is a nonlinear sigma model: 

= nig) - ^G^^^d.qd^q - ^G^^^d.rd'^r - ^(r, q) , (6) 



where TZ is the Ricci scalar of the 5d metric g^j,,^, r, q are scalar fields parametrizing the vector 
multiplets and hypermultiplets respectively. The metrics G'-'^-* and G^'"^ are given by: 

(1 -2g2)2 (2-3r2)2 ^ ^ 



and the potential V{r, q) is found to be: 

V{r,q) = (487r)2/3- 



2g2)(l-3/2r2)3/2 



The field r is one of the Kahler moduli, and can be thought of (at least when it takes positive 
values) as measuring the size of some 2-cycle in the compact space /C ^ The point r = 
therefore corresponds to a singular Calabi-Yau, and in this case the physics is still sensible 
because of the M2-brane wrapping states which become massless. 

Introducing the geodesic coordinates 

^ " ^y^arctanh (^ ^^2^2)3/2 / ' ^ " arctanh(^2g), (9) 

the action for the 5d effective theory becomes 

^5 = ^ y d'xV^nig) - J d'xV^ (^ld,<pd'^<P + \d,xd'x - V^(0, x)) (10) 



^When r is negative, it can be thought of as the volume of a 2-cycle in a manifold related to /C by a flop 
transition. 



6 



where = kI^^V^y = SttGs with Vcy the volume of the extra dimensions, G5 the 5D Newton 
constant, the 5d Einstein frame metric is g^^, and the effective potential is given by 

X) = ~g'^4>^')d + higher order (11) 

with the coupling = |(487r)^/^V^y ^ 1 in Planck units, since Vcy ^ in order that the 
SUGRA description is valid. 

In [38], the model (p-[8]) was used to demonstrate the possibility of conifold transitions 
occurring dynamicallyo Such a transition amounts to moving between the two vacuum 
branches: r = and g = 0,^ An important insight gained from the numerical studies of 
[38] is that the qualitative behavior of the fields and the scale factor(s) doesn't depend on 
the precise form of the potential. In particular it was insensitive to terms higher order than 
02^2 therefore expect that if we stay close to the origin, it is sufficient to study the 

minimal potential with no higher order terms: 

mx) = ^/0V (12) 

The property of ffTOl) of interest here is that the mass of one field depends on the ex- 
pectation value of another. This feature is not specific to conifold transitions, but arises 
whenever new massless states appear in the spectrum. For instance, if the Heterotic string 
is compactified on a torus at self-dual radius, a U{1) x f/(l) subgroup of the gauge group is 
enhanced to SU{2) x SU{2). Another example is a pair of closely separated D-branes. The 
ground state of a string stretched between them has a mass proportional to the separation, 
so as the branes coincide, unexcited strings starting and ending on different branes become 
massless. A similar phenomenon occurs in the small instanton phase transition [40|, and as 
indicated in |41| there are many other examples. 

Though the points in moduli space where (|T2l) applies are numerous, they are considerably 
fewer than the number of distinct vacua in the 'landscape'. So effects that trap the moduli 
near points like = x = may help to ameliorate the difficulties associated with extracting 
predictions from string theory. We now examine some of the special dynamics exhibited by 
motion in the potential (fT2|) . keeping in mind the problem of understanding how the fields 



^Dynamics near conifold transition points was also investigated in [52 . 

^"^The existence of two flat directions is not a generic property of singular Calabi-Yaus. Such a situation 
arises only when there are nontrivial homology relations among the collapsed cycles [49], A more typical 
potential contains terms oc 0"*, not divisible by 

^^Here we mean terms higher order in both and Xj not for example 4>^x'^ oi' X^- 
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and X acquire masses consistent with observations. 



3 Classical Moduli Trapping 

In this section we will think of ( ITOl) as a classical scalar field theory, minimally coupled to 
gravity. The dynamics is governed by a potential and by Hubble friction due to the expansion 
of the universe 

As can be seen from Figure [1], the potential restricts the moduli to a region of their 
parameter space consisting of two 'arms' (0 = and x = 0), and a central 'stadium']^ 




-1 -0.5 0.5 1 



Figure 1: Two views of the potential 1^(0, x) = |fi'^0^X^- The vacuum branches ('arms') 
= and X = 0, are thehorizontal and vertical axes. Following [36], we call the central 
region the 'stadium'. 

The effect of Hubble friction is that in an expanding universe, the scalar fields gradually 
lose energy, so over time their motion is constrained to smaller and smaller values of 0^x^- 
A natural question, then, is whether these two trapping effects, the potential and Hubble 
friction, are sufficient to mitigate the phenomenological obstacles mentioned in the introduc- 

Hubble friction is a somewhat misleading term. The force it refers to is proportional to the velocity of 
the scalar field, so a better mechanical analogue is air resistance. The distinction is important, since the 
fields cannot come to rest at a place where the slope of the potential is nonzero, as one might assume from 
the analogy with friction. 

^•^The terminology comes from [36] . 
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tion? This question was studied in [39]; here we extend the analysis taking particular care 
of the following two points. 

First, a crucial issue is that there is an ambiguity in many discussions as to what is 
meant by 'trapping' and 'stabilization'. It is sometimes assumed that all that is required to 
be consistent with observations is for the motion of the moduli to slow to a stop. In fact, 
one must also ensure that once they have come to rest, their couplings to matter are very 
weak, or they have acquired large masses. It is the latter requirement we pursue here. In 
other words, the eigenvalues of the Hessian of the potential must all be strictly positive (and 
sufficiently large). For a potential V = (f^x^ ■, we have: 

Since det = — 120^%^ vanishes for vacuum configurations, there is always one eigenvalue 
identically equal to zero, and so the phenomenological constraints cannot be satisfied. This 
is easily seen intuitively from Figure [H At any point on either vacuum branch there is a 
fiat direction the field can move in. The best one can do it seems, is to have Hubble friction 
keep the moduli some way down one of the arms, where only one of them has a mass. We 
will return to the viability of this option in section 13.31 

Second, we argue that the classical motion on the moduli space exhibits chaotic dynamics. 
As we will see, this implies that the classical stabilization found in previous works is only a 
pseudo-stabilization that traps fields for a finite duration, after which they leave the trapping 
region and continue their motion. 



3.1 Motion Without Hubble Resistance 

Before taking into account the effects of a cosmological background, it is instructive to study 
the dynamics of fIlOp in Minkowski space. In this circumstance X = 0; ^ = const, is a 
solution of the equations of motion, corresponding to constant velocity motion down one of 
the vacuum branches. And clearly no matter how long one waits, (f) will not return to the 
'stadium' area near the origin. One might guess that the same is true of almost straight 
trajectories where the fields move down one of the arms, but not quite along the minimum 
of the potential. However, it turns out that any momentum in the x direction, no matter 
how small, is sufficient to return to the origin. This was argued in [36] as follows. The 
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equations of motion are: 



4> + {9\')4> = (14) 
X + V)X = (15) 

where overdots denote derivatives with respect to time. If (p is large, and not moving very 
fast, then x is approximately a harmonic oscillator with frequency g(j). Following [36] we 
split the energy into two components, one corresponding to motion down the arm, the other 
to transverse oscillations: 

= = + \9^<P\^ (16) 

For any harmonic oscillator one can show that (x^) = E^^/g'^cfP', and replacing by its 
expectation value in (fT4|) results in an effective potential: 

0+^ = 0, Kfr = i?xM0) (17) 

For any positive E^, \4ff increases monotonically, thus returning to the origin. 



3.2 Chaos 

It has previously been noticed that the potential of Figure [1] has some of the features asso- 
ciated with chaos [HH]- For example, slightly different initial conditions can result in motion 
down different arms, and the effective logarithmic potential ensures that generically there is 
crossing of trajectories. Also since the potential in the stadium region is negatively curved, 
we can intuitively think of the resulting dynamics as a Newtonian approximation to geodesic 
motion in a negatively curved space, which is chaotic by the theorem of Horne and Moore 
[M] . The importance of these observations for the case of motion with Hubble friction war- 
rants a more concrete analysis of the appearance of chaoo , which we now perform using the 
Poincare 'surface of sections' method. The result is that indeed we find chaotic scattering 
in the stadium with almost regular motion in the arms. 



Chaotic behavior in cosmologies with scalar fields that have similar potentials was studied in [55l [56 
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Poincare Surface of Sections — Visualizing Chaos 

The Poincare surface of sections (henceforth PSS) technique was developed by Birkhoff [57] 
and Poincare and has been widely used in studies of chaotic dynamics. We will briefly 
review the idea of a PSS and then apply the technique to the problem at hand: motion in 
the potential V{(f),x) = \g'^4>^X^- 

The PSS is based on the idea that to observe chaotic behavior one need not examine 
the motion through phase space in all its complexity. Instead one can pick a codimension 1 
hyperplane in the phase space and mark each point that's crossed by a phase space trajectory. 
See Figure [2] for a schematic representation of the PSS technique. 




1 '1 I 



— ' ---t — ' 

Figure 2: Schematic diagram of a Poincare surface of section. The circles show the points 
where the phase trajectory cuts the hyperplane. Since the hyperplane is oriented, one can 
distinguish such points by the direction of the phase trajectory when it passes through. A 
PSS is a collection of such points with the phase trajectory moving in a specified direction. 
So this diagram really shows part of two PSS's — corresponding to the solid and hollow 
circles. 

It is a well known result of Kolmogorov [59], Arnold [60j and Moser [61] known as the 
KAM theorem that if a Hamiltonian system is integrable, then motion is confined to a torus 
in phase space. If the motion is bounded, then after a long time (many cuts through the 
hyperplane), a PSS traces out a planar section of a torus. In our example, the phase space 
is coordinatized by 0, XiP4>-,Pxi ^"^^ therefore 4 dimensional. By working on a fixed energy 
shell we reduce to 3 dimensions by solving for in terms of the other 3 coordinates and 
the total energy E. A codimension 1 hyperplane is then a 2d space, and if the system is 
integrable then the PSS is a smooth curve, i.e. a 'toric section'. 

We choose the hyperplane to be specified by the condition x = 0, and then plot the 
coordinates each time the trajectory pierces the x = slice in the > direction. In this 
way we build up a plot of successive {4>,p^) values. A smooth curve indicates the absence 
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of chaos. If on the other hand the system is not integrable, then we will see an irregular 
pattern of points. 

To develop an insight for the onset of chaos we consider the following potential which is 
related to (1121) by the addition of a control term ct(x^ + y^), 



For a ^ |(7 the motion is manifestly not chaotic as the potential is separable. In fact 
it turns out that cr ~ (7^ is sufficient to remove any indication of chaos from the PSS. To 
illustrate this, we will start with a = \g'^ = 10 and dial a down relative to g"^ so that we can 
see the signature of chaos emerge. We can see this effect clearly by comparing Figure [3] to 
Figure m where a has been reduced by a factor of 10. 



Figure 3: Surface of Sections for a = 10, l/2g'^ = 10, and E = 20. The points trace out a 
smooth curve. 



One should note that the smooth curve typical of an integrable system need not be an 
ellipse as in Figure [31 Figure [5] illustrates alternatives. The case of interest to us where a is 
exactly zero is more intricate because there are exactly flat directions in the potential. When 
the motion is far down one of the arms, we have seen that it is governed approximately by a 
harmonic oscillator in one direction, and a logarithmic potential in the other. We therefore 
see orbits that are approximately regular for a while, before degenerating into chaotic motion, 
as seen in Figure [H The almost regular orbits on the outer shell of the "onion" in Figures [6] 
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Figure 4: Surface of Sections for a = 1, 1/2(7^ = 10, and E = 20. Here the points fill out a 
region of the (j),p^ plane, indicating chaotic dynamics. 

represent the motion up an arm, which is always very close to a regular orbit. The irregular 
pattern of points filling up the central region indicates the chaotic nature of the motion in the 
stadium. Thus, though we may find stabilisation in the stadium using further techniques, 
we cannot use initial conditions to constrain exactly where in the stadium this stabilisation 
will occur. 

The significance of the alternation of chaotic motion with approximately regular orbits 
is that the trapping mechanism provided just by the potential is actually not effective in 
stabilizing or trapping any of the moduli. There might be a period of time where the fields 
are down one of the arms (so that one of them has a mass), but they will always return to 
the central stadium region before shooting off down another perhaps different arm. 

Next we will look at how Hubble friction changes this situation. Since the energy is no 
longer conserved once we turn on Hubble friction, a PSS analysis is no longer appropriate, 
however what we've learned in the above section will prove useful in all that follows. Specif- 
ically we will find that since the motion in the stadium is chaotic, even once the fields have 
lost a significant amount of energy to Hubble friction, they will always eventually scatter 
up an arm and it would thus require extremely special initial conditions to be trapped on a 
regular orbit. 
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Figure 5: Surface of Sections for a = 10, l/2g'^ = 10, and E = 20.0, for different initial 
conditions. 



3.3 Motion Including Hubble Resistance 

In an expanding universe, fields gradually lose energy. One might hope that consequently 
the fields eventually stop at some point away from the origin, thus leaving one field with a 
mass. We take the metric to be isotropic in 3 of the space dimensions and homogeneous in 
all 4: 

ds^ = -dt^ + a{t)Mx^ + b{tydf (19) 

As noted in the introduction, we choose this spacetime so as to compare our results directly 
with those of [381 EH ES]- Taking the gravitational background into account, the equations 
of motion become: 

+ {3Ha + Hb)^ + = (20) 

X+{3Ha + Hh)x + 9'<l>\ = (21) 

Ha + Hi - HaH, + ^kIT = (22) 
Hb + Hi - 2Hl + HaH, + '^KjT = (23) 

where T = ^ + j is the kinetic energy density of the fields. After some time, numerical 
simulations like those shown in Figure [7] suggest that the spacetime approaches a power law 
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solution: 

ds^ = -dt^ + e^dx" + e^^dy"" (24) 



for some constants q and q^. The energy density in the fields then satisfies: 

p=3/f„' + 3if.ff. = 5«(l±*) . (25) 

We expect that averaged over a sufficient interval, the kinetic energy is some fixed fraction 
of the total, though not necessarily 1/2 as for a harmonic oscillator. The velocity of the fields 
is therefore ~ ~ t^^, and the distance they have travelled ~ log(t). It follows that the 
fields move an infinite distance, and do not gradually approach any particular point in the 
moduli space. Moreover, since the scattering in the stadium is chaotic, after some time the 
fields will always eventually end up moving almost exactly down an arm, so in this sense 
Hubble friction does not trap the fields 

On the other hand, the effect of Hubble friction is to make the occasional sojourns down 
the arms rarer and longer lasting. By the time of nucleosynthesis (t ~ 10^^ Planck times), 
a typical journey down an arm may well last an extremely long time, so one might hope 

^^Were the dynamics not chaotic, one can imagine the motion confined to a stable orbit in the stadium 
which converges to the origin as the energy decreases. 



15 



_Q ^ I I I I I I I I I I I I I I I I I I I I 

10 100 1000 10000 100000 1e+06 

Time 

Figure 7: The scale factors as a function of time. After t ~ 1000 Planck times the spacetime 
is well approximated by power law expansion. Initial conditions here are (0, 0, x, x, a, a, b) = 
(0, 0.2, 0.2, 0, 1, 0.1, 1) and b is determined by the Friedmann equation. 

that over the time-scales we are interested in, the fields are essentially still. However, big 
bang nucleosynthesis (in particular the observed primordial Helium abundance) relies on a 
variety of processes that take place roughly between t ~ 1 sec, and t ~ 3 mins, and depend 
sensitively on the number of massless species of particle. The velocity of the moduli goes 
like 1/t, so with an initial velocity Vq, the distance in field space the moduli traverse during 
BBN is ~ folog(180) ~ S^o in Planck units. Since the mass of x say is gcp, we expect that 
over this time large changes to the masses of x aiid (p are possible despite the slowing effect 
of Hubble friction. Put differently, for the change in mass to be less than a GeV, the initial 
velocity must be fine tuned to Vq < 10^^^. The effect of Hubble friction is therefore still 
somewhat unpredictable^^ 

-"^^Even if the fields were effectively still over the course of BBN, there is no particular reason to believe 
that at the relevant time they should be in any particular place, i.e. near the origin, or down an arm. 
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Figure 8: x vs. (f) until t = 500 Planck times. Initial conditions here are {(j), (j), x, X, o, o, b) = 
(0, —0.2, 0.65, 0, 1, 1.1, 1) and b is determined by the Friedmann equation. The dotted line is 
the equipotential corresponding to the initial energy of the fields. 

It is worth noting that simulations performed in [39] up to t ~ lOOOtpianck suggest that 
the fields eventually get trapped near the origin. However, running the simulations for longer 
interval reveals that the apparent trapping at the origin can be misleading. Figure [S] shows 
the evolution of and x up to t = 500tpianck, when it appears that they are becoming trapped. 
Figure [9] shows the subsequent motion until t = 15000tpianck- During the intervening period 
the apparently stabilized field took a journey down the x = arm. Figure (TU] shows a 
different example of the same phenomenon, with the fields plotted against time. We can 
understand this as follows: with no Hubble friction, the fields are free to explore the full 
phase space between the lines of equipotential determined by the energy of the system. Once 
we turn Hubble friction on, the energy in the fields is no longer conserved and the available 
phase space decreases with time. However, because the arms are exactly fiat it is always 
possible for the fields to be sent arbitrarily far down an arm. Since the scattering is chaotic 
this eventually happens. The result is that instead of slowing the chaotic motion to a halt, 
Hubble friction merely stretches the motion out onto longer and longer timescales. 

These conclusions change significantly if one allows corrections to the potential. If 
the corrections do not lift the vacua (i.e. the additional terms are of the form 0™%" with 
m and n both ^ 0) then we still expect the fields to move down the arms every now and 
again, but there is no guarantee that they will return to the stadium. The derivation of the 
restoring force from the effective potential (|T71) depends crucially on the quadratic nature of 



17 



-0.6 -0.4 -0.2 0.2 0.4 0.6 

Phi 



Figure 9: x vs. (j) until t = 15000 Planck times. Initial conditions are as for Figure [HI 



If, on the other hand, a vacuum branch is lifted, say by the addition of a (f>^ ternj^, then 
Hubble friction indeed does have a trapping effect. Reduction in energy gradually confines 
to a smaller and smaller region of field space, and as t oo the energy vanishes, so there 
is some definite value of that the field configuration approaches. 



4 Quantum Particle Production 

Beyond the classical dynamics discussed in the previous section, essential quantum mechan- 
ical phenomena occur at ESPs. In particular, when a field is massless, on-shell excitations 
are produced at arbitrarily small energy cost. 

Particle production at ESPs has been studied previously in the context of preheating 
models [^21 E31 EH ESI ES] , and also in various string theory scenarios [571 ESI SD SH ESI 
[701 [7n [12]. But it is of interest here because it results in a potential which confines the 
fields near the extra species point [41j. The rough idea is that as one of the fields becomes 
massless, a burst of on-shell particles is produced. The mass of these particle increases if the 
fields subsequently move away from the ESP, so energetics force the fields to stay nearby!^ 

^^This is precisely what one expects in the neighborhood of an arbitrary conifold point, i.e. one that is 
not necessarily a transition point between two topologically different families of Calabi-Yaus. 

^^Particle production is not the only quantum modification to the classical motion. There are also correc- 
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Figure 10: The fields and x &s functions of time. Looking only up to t ~ 10^ one migfit 
tliink tlie fields were trapped at tlie origin. But tlien at t ~ 2 x 10^ starts to increase. Initial 
conditions used here are: (0, 0, x, Xi c-; 0,1 b) = (0, 0.2, 0.2, 0, 1, 0.3, 1) and b is determined by 
the Friedmann equation. 

For example, consider the potential (fT2l) . and for initial conditions take and x approx- 
imately homogeneous in space, positioned along one of the vacuum branches, and moving 
towards the origin. For definiteness say x = 0, > 0, and < 0. The effective mass of 

tions to the kinetic terms and loop corrections to the potential (Coleman- Weinberg corrections). There is 
however a set of circumstances in which particle production is the most important effect. It was shown in 
[44] that for weak coupling {g^ ^ 1) the effect of on-shell production dominates over propagator corrections 
from virtual scattering. This will be the case here as well, since the coupling g must be small to ensure the 
validity of the SUGRA description. The Coleman- Weinberg potential requires a bit more care, since after 
SUSY breaking loop corrections should be expected to play an important role. In [41] it was argued that 
these corrections are subdominant because on-shell production always gives a positive contribution for both 
bosons and fermions, whereas fermions running in loops will give an opposite contribution to the effective 
potential from that of bosons. If SUSY breaking is soft, this implies that the induced potential from on-shell 
production should dominate the effective potential relative to the contribution from loops. However, a better 
understanding of the effects of SUSY breaking remains an important goal for string cosmology. 



19 



the X particles is = which vanishes as crosses the origin. So a burst of on- 

shell X particles is produced. Now can no longer move freely away from the origin since 
conservation of energy places a restriction on how massive the incoherent x excitations can 
become. This results in a potential, trapping cj) at the origin, which gets steeper as more x 
quanta are produced. Since there is a interaction term, scattering of x particles results 
in excitations of the (f) field as well, effectively driving both fields towards the origin. 

The aim of this section is to study how particle production might change the conclusions 
of the classical analysis of the previous sections. We use the same metric ansatz as before: 

ds^ = -df + a{ty {dxj + dxl + dxl) + h{tfdxl (26) 

with stress tensor Tj} = diag(— p, p, p, p, p). And for simplicity we take the fields to be moving 
along the x = ^irm initially. The equations of motion are then: 

^+{3Ha + H,)^ = -V'{<j)) (27) 
Ha + Ha i3Ha + H,) = '^{p - p) (28) 
H, + H, i3Ha + H,) = ^{p- p) (29) 

where Ha = d/a and H^ = b/b are the Hubble parameters. To study excitations of we 
substitute x = x{t) + '^x(^? u) iiito the x equation of motion ([21]), which results in: 

5x + {3Ha + H,) 5x - (^^i + -^^l - 9'<Pit?^ = 0, (30) 
A slightly modified Fourier transform gives a simple mode equation: 

6x= [ d'^k dh^^e"'-'^+"''y (31) 
J V a% 

Xk + ujl{t)xk = ^ (32) 
where the frequency ujk of the k^^ mode is: 

'^'W = ^ + § + 9'^\t)-\{9Hl + Hl + QHaH, + QH^ + 2Hy) . (33) 



To proceed we need to compute the contributions to the right hand sides of equations (I2714291) 
(i.e. V\ p and p) resulting from a burst of particle production. This can be determined from 
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n^, the number density of particles produced with physical momentum k, in the following 
way: 

^^ n^ = number density in all modes (34) 

^^--^u;«;n« (t^K = uJk with {k, /C4) = (a/?, 6^4)) (35) 

To find the potential it is useful to recall the case with 2 fields: V = |5'^</>^X^! and approxi- 
mate ~^ (x^)- For the incoherent x oscillations associated with production of x particles, 
we have: 

g^(j)^ g\(f)\ 

since for slowly varying (j), x is approximately a harmonic oscillator with frequency g(l). In 
the last step we used ~ gn^\(j)\, which amounts to approximating the frequency by the 
contribution to it from the effective mass onlyl^ 
then given by: 



The derivative of the potential V'{(j)) is 



V'{<t>) = 9' ix') = gn^^^ (37) 

Finally, the pressure p can be found by solving p = — {SHa + Hb) {p + p), and using the 
equation of motion. 

It is important to note that even if there is just a single production event, is a function 
of time, because the particles are diluted as the universe expands. Since ra^ is appreciable 
only for A; <^ 1 in the regime where effective field theory is useful, the produced excitations 
behave like dust: 

where t# is the time of production. The number density per mode ni^{t^) was calculated for 
an isotropic background in |41j, and the results generalize (as derived in the appendix) to 



^^In other words, we assume ~ g'^cj)^ + f{k), and take f{k) to be small. This is a good approximation 
because f{k) is suppressed by powers of time, and we expect the time of production t# to he ^ tp because 
natural initial field values are a distance ^ 1 from the origin, but initial velocities must be ^ 1 to justify an 
effective field theory description. 
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our 5d follows 



1 • 1 • 

P = o<f + Px= o4>'^i^)+ 9n^it#) 



(39) 



a\t)b{t) 



(40) 



(41) 



(42) 



where v# = -i\t=t^. 

4.1 Numerical Simulations and Comparison with Classical Effects 

Our aim now is to understand how production of on-shell excitations changes the results of 
section I3.3[ Rather than work with typical initial conditions, we consider the 'worst case 
scenario' from the point of view of trapping, i.e. motion directly down one of the arms of the 
potential (classically there is no trapping at all: 0(t) ~ log(t)). Through this analysis, we 
intend to set a lower bound on the efficiency of trapping. With this in mind we include only 
a single particle production event, and do not alter the potential in subsequent crossings of 
the origin. 

More general initial conditions do not allow the fields to move arbitrarily far from the 
origin. And it is argued in [41j that successive particle production events are enhanced by 
parametric resonance since they take place in a bath of the particles already produced. 

It is not easy to guess the net result of including particle production, since there are 
several competing effects: 

• Hubble friction slows down the motion of the fields 

• The expanding background dilutes the produced particles thus reducing the slope of 
the induced potential. 

• The particles backreact on the background 

In jH] the first two of these effects were taken into account, but the scale factor was 
assumed to be a power of time. Here we include all three effects. For typical initial conditions 



'In D + 1 dimensions: oc dkk^ exp (~7-|~ ) ° , where 7 is a number of order 1 
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Figure 11: Evolution of with time for typical initial conditions. 



(with either a 4d or a 5d background), we find solutions like those in Figure [TTl A plot 
of log(0) against log(t) (Figure [12]) reveals a power law decrease in the amplitude of 
oscillations, suggesting that after is restricted to an ever-decreasing range of values, it 
approaches like a negative power of t as t oo. Since this is a lower bound on trapping 
efficiency, we expect that motion with generic initial conditions (not along an arm of the 
potential) has the same features. This is to be contrasted with the classical situation of 
Figure [TDl where the fields are not contained in an envelope with monotonically decreasing 
amplitude. Quantum particle production thus provides a viable mechanism for trapping 
some of the moduli at points where extra massless species appear. 



5 Conclusion 

We have considered the effects of additional light degrees of freedom on both the classical 
and semi-classical dynamics. In the purely classical case we have illustrated the appearance 
of chaotic scattering in the stadium region, and found that, somewhat surprisingly, vestiges 
of the chaotic dynamics are important even when Hubble resistance is taken into account. 
In particular, rather than slowing the fields to a halt, the effect of the expansion seems 
to be to stretch the motion out on ever larger timescales. The point is that the chaotic 
dynamics ensures that the weakest link in the stabilization mechanism will, sooner or later. 
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Figure 12: Here the maxima of are plotted on a log-log scale, showing the power law 
decrease of the amplitude. 

be probed. And because the potential contains flat directions, these loci form the weak 
links and they allow fields to roll away from any would-be stabilization location. Moduli 
stablization therefore seems unlikely to result from classical motion in the quartic 0^%^ 
potential. Higher order corrections may change this conclusion significantly, but the end 
result (i.e. where the fields finally rest) will still depend sensitively on initial conditions if 
the corrections do not lift the flat directions. 

These caveats about stabilization do not apply once the effects of quantum particle 
production are taken into account. The amplitude of the expectation values of the fields 
decreases like a power of time. Though we only simulated the effect for the field, once 
incoherent x excitations are present, the term allows scattering to produce particles 
also, which produces a potential for x similar to that for 0. 

An important problem that remains is to consider additional quantum effects. In this 
semi-classical analysis we have only focused on on-shell particle production. However, un- 
derstanding propagator and loop corrections in the absence of supersymmetry is a vital next 
step. 

One obstacle to using the above framework to stabilize moduli is that in the known 
models whose low energy theory includes a piece of the form ffTOj) . only a proper subset of 
the moduli possess interaction terms of the type required. For example in a Calabi-Yau 
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compactification of M-theory, if we approach the conifold point by blowing down 2-spheres, 
then in general only the Kahler moduli will couple to the brane wrapping modes. Separate 
considerations, such as the inclusion of p-form flux windings or nonperturbative effects like 
gaugino condensation, are still needed to produce a phenomenologically acceptable theory. 

One may also wonder whether the linear potential resulting from particle production 
really solves the phenomenological problems mentioned in the introduction, since it is not 
a mass term. However, the concept of particle number is only well defined far away from 
the time when crosses the origin, so it seems likely that the sharp cusp in the potential is 
smoothed out. A large mass term is then present, but large interactions also. 

Finally, it is important to note that the potential arising from particle production is 
proportional to the number density of the excitations. So as the universe expands, the 
trapping effect becomes weaker and weaker. Therefore, while it may help to ameliorate 
the issues moduli create with primordial nucleosynthesis, late-time effects such as varying 
couplings are still a problem. 

The important case of particle production and dynamics in a contracting universe has not 
been previously studied and will be the subject of a future paper [72j. It is not immediately 
obvious what the effect of contraction will be and it may be of great relevance to models of 
the universe that include a phase of contraction such as the cyclic universe |73j . 
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A Particle Production 



A formal solution to the mode equation ( l32l) is given by 



^^^t) = ^^^^^ .-^r-.jt'W ^ Pkjt) ^.p^^^ew (43) 

Normalizing Xk by using the Klein-Gordon inner product is then equivalent to |a/c()f:)p — 
|/3fc(t)P = 1, which can be used to write the equation of motion as 

Pk = ^e-'^f'-^^'^'^'a, (44) 

If initially we have <C 1 and ~ 1 then at late times, the number density of particles 
produced in the fc*^ mode is Uk = \Pk\^- Assuming that Pk is always small, the production 
can be estimated by: 

(3,^ f dt^e-''^'-^^''^''' (45) 

J ^Uk 

We assumed earlier that the particle production happens in a short burst as (j) crosses the 
origin. To justify this approximation, note that the integrand of fHSl) is only significant when 
Uk ujk, or in terms of dimensionless quantities when: 

^i>l . (46) 

To get an idea of how ^ behaves, consider the following analytic solution of the equations 
of motion 

4>{t) = v#t#\og(J-^ (47) 
Ha{t) = f (48) 
Mt) = (49) 
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Figure 13: The non-adiabatic parameter ujjij^ for very long wavelength modes (plotted here 
for k = 0). The other parameters are t# = 10^ and gv^ = 10~^^, consistent with the effective 
field theory approximation. The sharpness of the peak increases with and is independent 
of the choice of Kasner parameter q. 



which results in the following mode frequencies: 



^2 



t2q t2(l-3q) 4^2 



+ 1^ + a'vltl 



log — 



(50) 



A natural initial field value is ~ 1, while ^ 1 in order that an effective field theory 
description is appropriate. It follows that the time t# when = is large. Figure dSl shows 
a plot oi Lo/uj'^ for typical values of and from which one can see that indeed a burst 
of particle production does occur at t = Less obvious is that lo/uo'^ is large near t = 
as well. This particle production is due to the expansion of the spacetime (essentially the 
l/4t^ term in (1501) ) rather than time dependence of 0. Here we will only calculate the effect 
of the latter, i.e. production at t ~ 

For a general solution, the mode frequencies are: 
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Figure 14: The deformed contour for evaluating f3k- 
integral we evaluate. 



The arc near is the only part of the 



and neglecting the gravitational terms ( since ^ ^pianck) gives; 

We now estimate the integral fHSj) following the method of \7^. Note first that the 
integral is dominated by times when uj is small. By allowing t to take complex values, one 
can take advantage of the analytic structure of Uk to compute the integral. In particular, 
one can deform the contour of integration so as to pass near a zero of ook, located at a 
(necessarily complex) time t# whose real part is close to t#. Since Uk is real for real values 
of t, the zeros come in pairs related by complex conjugation, and from the single-valuedness 
of ul, we see that t# and t# are the branch points of a square-root branch cut, which for 
convenience we take to pass through too rather than along the real axis. The idea is that we 
will approximate jSk by integrating around an small arc near that interpolates between 
the two steepest descent paths, as in Figure [HI 

On this path we can expand Uk in a Taylor series around t^: whose terms are half-integral 
powers because of the branch point. In general, if f{x) is analytic in a neighborhood of x = 
and /(O) = 0, then we have: 



X 



r ' ' 4 V3 ' ' 8 /'(O) 



(52) 
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The frequency therefore becomes: 



(53) 



and integrating with respect to time: 

ujk{t')dt' = I 

^ 

Substituting the above expressions into equation (H5ll . one finds 



fe = e 



2i !l* u^^(t')dt' (If d5_=^,/J{u)s^/^ 



(54) 



(55) 



where 6 = t — the C denotes the red arc in Figure [T3] interpolating between the two 
steepest descent contours. The bracketed integral can be performed by a change of variables 



^^ = 1^ V'/(^*)'^^^^ which (see [71]) closes the contour C to a loop as well as simplifying the 
exponential. The result is: 



27r 



(3k = ^^P ( / ujk{t')dt' ) exp I —2i j uJk{t')dt' 



(56) 



The first integral is real, contributing only a phase to the exponential which will cancel in 
l/^fcp. The second integral can be approximated j^^ ujk{t')dt' ~ yIm(t=K)c<jfc(t#), where 7 is a 
constant of order 1. The result is: 



TT 



37lm(t*)u)fc{t#) 



(57) 



Thus, to calculate the leading contribution to particle production it is sufficient to identify 
the real and imaginary contributions to the zeros of uj^ Hll For small /c, we expect — to 
be small, so the following approximation is useful: 



= U^liU) ^ Co + Ci{U - t#) + C2{U - t#f 



(58) 



^"'^We only get a sensible answer by deforming the contour into the lower half plane. This can be shown 
[74] to be a consequence of choosing the positive square root in ((53|l . 
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where the Taylor coefficients are: 



Co 



Cl 



^4 



C2 = 0{k')+gY{t#) 



(59) 

(60) 
(61) 



Here we have used 0(t#) = 0. We can now find the zeros of Uk by solving (1581) . As argued 
in more detail in [71], we are interested in the zero in the lower half plane: 



-Cl - a/ci - 4coC2 

2C2 



(62) 



For small k (so that terms 0{k^) can be neglected), the comoving number density per mode 
is: 



^'^-g7lm(t.Vfc(i#) 



71 



-7Co/v^ 



TT 



(63) 
(64) 



The total number density (in physical rather than comoving coordinates) in all modes is 
then found by integrating: 



1 



a3(t#)&(t#) \3J J i2n 
14472 



^ r d^k ( 1 
exp — 7- 



a2(t#) 62(i^) 



(65) 
(66) 
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